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Abstract 

We obtain an exact finite-size expression for tlie probability tliat a percolation hull will touch 
the boundary, on a strip of finite width. Our calculation is based on the g-deformed Knizhnik- 
Zamolodchikov approach, and the results are expressed in terms of symplectic characters. In the 
large size limit, we recover the scaling behaviour predicted by Schramm's left-passage formula. We 
also derive a general relation between the left-passage probability in the Fortuin-Kasteleyn cluster 
model and the magnetisation profile in the open XXZ chain with diagonal, complex boundary terms. 

1 Introduction 

(N 

Percolation models in two dimensions play an important role both in theoretical physics and mathematics. 
On the physics side, it was one of the first models where the Coulomb gas approach [T^ was used 
to predict the critical exponents [23 , and where the concepts of boundary conformal field theory (CFT) 
fH were put in practice [4]. Nowadays, it still attracts the commmiity's attention, especially for its relation to 

logarithmic CFT. On the mathematics side, many rigorous studies of percolation have been pursued |28j, 

Sand Smirnov proved [26. that site percolation on the triangular lattice has a conformally invariant scaling 
I I limit, described by Schramm-Loewner evolution (SLE) with k = Q. Also, in combinatorics, the Razumov- 

Stroganov relation pi 12011^ identifies the components of the percolation transfer matrix eigenvector with 
the enumeration of plane partitions and alternating sign matrices. 
^ The main objects of study in percolation are the percolation clusters and the lattice curves surrounding 

them, known as hulls. In the scaling limit, the correlation functions of these hulls are conjectured [53] 
to be described by a Coulomb gas CFT [TH (TH] , and thus to satisfy some partial differential equations 
(PDEs) given by the "null-vector equations". Some of these PDEs can be solved explicitly, e.g., the one 
for the crossing probability (the probability that a cluster connects two sides of a rectangle) . A very 
fruitful approach to relate CFT and SLE is to express the null-state equations of CFT as martingale 
CS| conditions for the SLE observables [2 [S] . 

The left-passage probability Pic{t{z), i.e., the probability for an open, oriented hull to pass to the left 
^ of a fixed point z of the system, is one of these observables that can be easily obtained for percolation 

r ""j (and more generally, for the Potts and 0{n) models) both from the CFT and SLE viewpoints. In the SLE 

rS literature, this result is known as Schramm's formula 24 . In this paper, we address the determination 

of Picft on the lattice, in the infinite strip geometry, using rigorous techniques based on the Yang-Baxter 
and quantum Knizhnik-Zamolodchikov (gKZ) equations, as well as the Bethe Ansatz for the related 
six-vertex model. 

The (jKZ approach is particularly powerful for loop models with a trivial partition function Z — 
1 [II [TOl HU [121 [30] . In several cases, it allows the explicit determination of the dominant transfer matrix 
eigenvector, and it turns out that the components of this vector are integers enumerating plane partitions 
and alternating sign matrices. Also, this technique was recently used by one of the present authors 
to calculate a finite-size correlation function in percolation, namely the "transverse current" across a 
strip [B]. A complementary approach, with a larger scope, is to map a loop model onto an integrable spin 
chain 2J, and use the Bethe Ansatz to obtain correlation functions in the form of determinants 17, 18]. 
This approach actually extends to the Fortuin-Kasteleyn (FK) cluster model with cluster weight Q, and 
the percolation model is recovered for Q = 1. 



* Yacine. Ikhlef@unige.ch 
t Anita. Ponsaing@unige.ch 



1 



The layout of this paper is as follows. In Section [2] we recall the exact equivalence [2] between bond 
percolation on the square lattice and the Temperley-Lieb loop model with weight n = 1, briefly review 
its conjectured relation to SLEg, and state our main results. In Section [3j we set up our notations for the 
transfer matrix and recall the basic steps of the qKZ approach. In Section |4] we derive explicitly, for finite 
strip of width L, the probability Pieft(z) with z on the boundary of the strip. For a homogeneous system, 
we obtain fractional numbers with a simple combinatorial interpretation. In the large-L limit, we recover 
the power law predicted by CFT and SLE. In Section [5j we generalise to a generic point z: it turns out 
that similar symmetry and recursion relations hold, but are very difficult to solve in practice. However, 
we obtain two promising results in this case. First, we calculate Pie{t(z) numerically for homogeneous 
systems with up to L = 21 sites, and observe good convergence to Schramm's formula. Second, we prove 
that, in the FK model, -Pieft(2) relates very simply to the magnetisation profile in an open XXZ spin 
chain with diagonal, complex boundary terms. We give our conclusions and perspectives in Section [6j 



2 Percolation, Temperley— Lieb loops and SLE 
2.1 Percolation hulls and their scaling exponents 




Figure 1: Left: the original lattice C where bond percolation is defined (black) and its medial 
lattice Ai (grey). Right: an example cluster configuration (thick lines and dots) and the asso- 
ciated loop configuration (thin lines). On the left (resp. right) boundary, all edges are empty 
(resp. occupied). 



Consider a square lattice £, on which each edge can be occupied by a bond with probability p, or 
empty with probability (1 — p). The connected components of the graph formed by all the sites and 
the occupied edges are called percolation clusters. We now look at the medial lattice A4 formed by the 
mid-edges of £, where a loop configuration is associated to each cluster configuration ^2j (see Figure [T]). 
These loops follow the external boundaries and the internal cycles of the clusters, and are called the 
percolation hulls. The model describing these loops is called the Temperley-Lieb loop model. 

As usual, the scaling limit is defined by fixing a domain fl of the plane, and covering it by a square 
lattice with spacing a — 0. The scaling properties of percolation hulls at the critical point Pc = 1/2 can 
be determined by the Coulomb gas approach, yielding the £-\eg "watermelon" exponents Xg — (£^ — 1)/12, 
the fractal dimension df — 7/4, and the correlation-length exponent ly = 4/7. 



2.2 Strip geometry and the SLE model 

We suppose that the system is defined on an infinite strip of A^, of width L, where L is odd. Let us 
specify the boundary conditions (BCs) which shall be used throughout the paper. On one boundary, we 
set all the edges of £ to be occupied (wired BC), and on the other boundary, all edges are empty (free 
EC). This simply means the hulls must reflect on both boundaries. Since L is odd, there exists an infinite 
open hull propagating along the strip, which we shall denote 7 (see Figure |2]). 

In the scaling limit, the random curve 7 is conjectured to be distributed as chordal SLE^^g- In the 
SLE model, the left-passage probability Pieft is obtained by solving an ordinary differential equation of 
order two. If we normalise the width of the strip to L x a = 1 and denote hy x E [0, 1] the horizontal 
coordinate, Schramm's formula reads [24] 

1 r(4/K) ^ /I 43. 



Pioft(2:) = ^ t8-K\ cotan -Kx X 2F1 -, -; -; -cotan nx ] , (2.1) 

^ i (^j \Z K / 
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Figure 2: (a) A configuration contributing to the boundary passage probability P\j on the selected 
section (dotted line), (b) A configuration contributing to the boundary passage probability Ph- 
In both cases, the dots indicate for which edge the passage probability is defined. 



where 2-^1 is the hypergeometric function, which gives the probabihty of "touching" the boundary 

r(4/«:) 



-PiGft(a;) 



(2.2) 



.0 2V^r(^) 

with K = 6 for percolation. In the present work, we derive some exact resuhs for Picft(a;) in the lattice 



model, i.e., we look for analogs of (2.11 and (2.2) infinite size. 



2.3 The Fortuin— Kasteleyn cluster model 

Most of the above results can be generalised to 4 < k < 8 by considering a modified cluster model, called 
the Fortuin-Kasteleyn (FK) model, where each cluster gets a Boltzmann weight Q (the critical regime is 
< Q < 4) , so that the Boltzmann weight of a cluster configuration C is 

'W^[C'] — ^,#occupicd cdgos(C) 

Loops are defined on the medial lattice similarly to percolation, and, using the Euler relation, one finds 
that the above Boltzmann weight can be written as 



W{C\ oc 



T^/^closcd loops(C) 



)7t^occupicd cdgcs(C) 



This defines the Temperley-Lieb loop model with weight n = ^/Q (see Section 3.1 1. It is conjectured 
(and proved for Q = 2) that the hulls of FK clusters are distributed in the scaling limit as SLE^ with the 
relation 



~2 cos ■ 



An 



4 < K < 



2.4 Statement of results 

• In the percolation model with n = —{q + q^^), using the gKZ approach, we obtain explicitly (see 
Section 4.4.1 ) the probability that the open path 7 passes through a boundary edge. It reads: 



Pb{zi,Z2, ■■■,Zl) 

Pb{w;zi,Z2, ■■■,zl) 



XL-i(z|, ■■■,zl) XL+iizf, zl,zl,...,zl) 



,,2 ,,,-2^ 



[q ^-q){w'^~w ^xl+jW^I,...,zI) XL+i{{q/w)'^,zl,...,zl) 



Ut=iH^/w,z,) 



Xl{zI,...,zIY 



depending on the sublattice where the boundary edge sits (see Figure |2|. In the above expressions, 
q has been set to exp(2i7r/3), corresponding to n = Q = 1, the zfs are the vertical spectral 
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parameters, w is the horizontal spectral parameter, xl is the symplectic character (see notations 
in Section |3| , and we have defined 

k{a, b) ~ {q/af + [a/qf - 6^ _ 5-2 



For a homogeneous percolation system (see Section 4.4.21, this becomes 



Ph = 



Av{L) Av{L + 2) 
N^{L + 1Y 



and 



Ph = 



A{Lf 



4^ N^{L+IY Av{LY 



where A{L), Av{L), and N%{L) are the number oi Lx L alternating sign matrices, L x L vertically 
symmetric alternating sign matrices, and L x L x L cyclically symmetric self-complementary plane 
partitions respectively. Moreover, for a large system size L, both Pf, and Ph scale like L~^/'^ (see 



Section 4.4.3), which is consistent with (2.2 1 



• In the critical Fortuin-Kasteleyn model with generic parameter Q e [0, 4], we define the probability 
Xj that the path 7 passes through the jth horizontal edge of a given section as in Figure [2^ (so 
that Pi, ^ Xl). Defining Picft on the dual of A4 , one can write 



-Pioft {xj+1/2) - Picit {xj-1/2) = (-1)^ , where 



Xj := j/L. 



We find the relation (see Section 5.3 ) 



(*o|*o) 

where j^E'o) is the groundstate eigenvector of the open XXZ Hamiltonian 

L-l r 



where ^/Q = -{q + q^^). 



1 



{q~q-'){<jt~al). 



3 The gKZ approach 

3.1 The Temperley— Lieb loop model 



G 









Figure 3: The lattice with two reflecting boundaries. 

The Temperley-Lieb loop model with wired (or reflecting) boundaries [TTJ [T^l [131 [TS] is defined on a 
square lattice, where each face is decorated with loops in one of the following two ways: 



and 



Every closed loop gets a weight 



The chosen boundary conditions for this model require that the lattice is infinite in height and of finite 
width L. On the left and right are refiecting boundary conditions, as in Figure [3j 
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Drawing a horizontal line across the width of the lattice, we consider the connectivities of the loops 
below the line while ignoring the paths the loops take, as well as any closed loops. We refer to this 
pattern of connectivities as a link pattern, and we denote by LPl the set of link patterns for a given 
system width L. For odd system size L = 2m — 1 they are enumerated by the mth Catalan number, 
(2m)!/(m!(m + 1)!). An example link pattern for L = 7 is 



We label the link patterns by \a), using the shorthand notation of "(•••)" to indicate a pair of sites 
connected by a loop, and "|" to indicate the single unpaired loop which always exists in an odd-sized 
system. As an example, the above link pattern is indicated by \a) = |(())|()). 

The link patterns for a fixed L form a representation of the Temperley-Lieb algebra, generated by 
{ej, 1 < i < L — 1}, with e, depicted as 

(■ i+l 

The relations for the Temperley-Lieb algebra are 



n ej 



if \i - j\ > 1. 



(3.1) 



Multiplication corresponds to concatenating the depictions of the generators, giving a weight of n to 
every closed loop, and disregarding the paths the loops take. In this way, we obtain relations between 
the link patterns such as e4|(())|()) = ||()()()). 
A state in = span(LPi) is written as 



aeLPi, 



We look at all the possible configurations of two rows of the lattice, and consider how they send a given 
link pattern to another. We can write this as a matrix t which acts on V^, and we refer to this as the 
transfer matrix. 

We take an arbitrary initial state |in) and act A'' times with the transfer matrix t. As N ^ oo, we get 

lim |in) oc A^ I*) , 

AT— >-oo 

where A is the maximum eigenvalue of t and l^") is the corresponding eigenvector, also known as the 
ground state. When n = 1 all the weights are probabilities and therefore A = 1. The components of 
can be thought of as the relative probabilities of the possible link patterns. In a similar way we define 
which is the groundstate of the rotated lattice, giving the relative probabilities of upward link patterns. 
The inner product between upward and downward link patterns is simply 



a) 



#closed loops 



Va,/3. 



Hence, the expectation value of some observable O reads 

(^-lOI*) 



(*l*> 



3.2 The i?-matrix 

The possible states of each lattice square are described by the i?-matrix. 
Definition 3.1. 

_ [qz/w] 



R{w, z) = w ■ 



[qw/z 



+ 



[z/w] 
[qw/z^ 



where 



[z] ■.= z- l/z . 
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Lemma 3.1. This R-matrix satisfies the Yang-Baxter equation (YBE) 





(3.2) 



the unitarity relation 



and the crossing relation 




(3.3) 



w ^ qw 



(3.4) 



z z 

Definition 3.2. The corresponding operator, acting on Vj,, is 



k/w] ^ _ \w\_ 

[qw] [qw] 



(3.5) 



3.3 The transfer matrix, symmetries and recursions 

We now define the transfer matrix t, which describes aU the possible configurations of two lattice rows 

mm. 



Definition 3.3. 



t{w] 01, ... , zl) = Tru, [R{w, zi) . .. R{w, zl)R{zl, l/w) . . . R{zi, l/w)] 



or pictorially, 



t(w; Zl, . . . , Zl) 



} 



Lemma 3.2. Thanks to the YBE (3.2), the transfer matrix satisfies the interlacing relation 



Ri{zi/zi+i)t{w; Zl, Zi+i) = t{w; Zi+i, Zi)R^iz^/ z^+i), 



(3.6) 



pictorially, 



c 



) ■ 



Considering the possible configurations of the two tiles at either position 1 or position L of the transfer 
matrix also gives us, respectively, 



t{w;zi,Z2...) ==t{w;l/zi,Z2,...), 
t{w; . . . ,zl~i,zl) = t{w; . . . ,zl-i, 1/zl). 



(3.7) 



Lemma 3.3. By acting the transfer matrix on a small link from site i to i + 1 (denoted by ipi) and setting 
Zi+i — qzi we find the relation 



tL{zi,z,+i ^ qZi) o ipi = (PiOtL^2{h,'Zi+i), 
where z means that z is missing from the list of arguments. 



(3.8) 
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(3.31 



Proof. Considering first the bottom row, we use the crossing relation (3.4 1 and then the unitarity relation 



and see that the bottom row no longer depends on Zi and Zi+i. We repeat the procedure for the top row: 

Zi qZi Zi Zi 



w 

and the result follows. 



□ 



3.4 The dominant eigenvector 

It is possible to show (see for example fT^) that two copies of the transfer matrix with different values 
of the parameter w commute, and therefore that the groundstate eigenvector does not depend on w. 
Explicitly, the eigenvalue equation thus becomes 



t(w] zi, . . . , zl) \^{zi, zl)) = \^{zi, zl)) , 
with the ground state eigenvector given by 

. . . = ^ ll)a{zi,...,ZL)\a) . 



(3.9) 



From the expression of the i?-matrix, the coefficients in the eigenvalue equation ( 3.9 1 are all rational func- 
tions of the Zj's, and hence one can normalise |5'(zi, . . . , zl)) so that all the components 4'aizi, . . . , z^) 
are Laurent polynomials in the Zj's. Moreover, one requires that these components have no common 
factor. 



With this normalisation, the interlacing relations (3.6) and (3.7 1 yield the q-deformed Knizhnik- 



Zamolodchikov equation for the ground state eigenvector, expressed in the form 

R,{zi/zi+i) \'i'{zi, . . . = |^'(2:i, . . .,zl)) , 

|*(zi,...,zl)> - |*(1/zi,...,zl)) , 
|*(zi,...,zl)) = |*(zi,...,l/zi)), 

where 7ri/(z^, z^+i) = /(z^+i, z^). 



Acting on |'I'L-2(^i, ^i+i)) with both sides of (3.8), we get 



tL{zi,Zi+i = qZi) ipi \^L-2izt, Zi+i)) = ifii |*L-2(zi, z,;+i)) . 

Since the ground state is unique, |\I'i(zi_[.i = qzi)) and ^pi |'I'L-2(zi, Zi+i)) are linearly related, and one 
can show that \'^) satisfies the recursion relation 



with 



|*L(zi+i = Qz,)) = |^'L_2(zj,Zi+i)) X (-1)^ Yl k{z„Zj), 
k{a, b) = [qb/a] [q/ab] . 



(3.10) 



3.5 Solution for the eigenvector 

The gKZ equation forces certain symmetry requirements on the components of |^), which lead to the 
solution for the ground state. For instance, 

^i(...o...) =(-i)*(*+i) n fc(z„zo n 

l<«<i<^ ^<i<j<L 

By considering the Dyck path representation of the link patterns, one can write the other components in 
terms of factorised operators acting on this component. We will not give the explicit solution here as it 
is not needed for our calculations. A full explanation of the procedure is in Section 4.1 of [TT] . 



7 



3.6 The normalisation factor Zi^ 

The normalisation factor Z/, is defined as 

. . . , zl) := ^ ■0a(^i, 



To express Zl, we first define the completely symmetric polynomial character of the symplectic group. 
Definition 3.4. The symplectic character x\ associated to a partition A is given by 



where Sj = L — j + 1 and fij = Xj + Sj . 

Throughout this paper, we shall restrict to the partition Xj = 
recursion 



det 




det 





L-j 
2 



, for which x has the special 



n.=,., = (-1)^ n xr "(•■•, 

We will use the shorthand notation xl(- • •) '■= X\\- ■ Oi with the particular choice of A given above 
Proposition 3.4. The normalisation Zj^ is given by 

Zl = Xl{zI, . . .,zl). 



(3.11) 



Proof. The recursive property ( 3.10 ) of the ground state eigenvector is easily extended to its components, 
and thus to the normalisation. 



ZL{zi+i = qz.,) = ZL-2{zi,Zi+i) X (-1)-^ Y\. Kzi,zt) 



(3.12) 



As Z]^ is a symmetric function (easily proven using (3.3) and the qKZ equation), this can be generalised 
to 

ZL{zj = qzi) = ZL-2{zi,Zj) X (-1)-^ Y\_ k{zi,Zi). 

The symplectic character xl{zi, ■ ■ ■ , z\) also satisfies these recursions, and it is straightforward to show 
that these recursions are enough to satisfy the degree of , which is set by solving the gKZ equation. It 
remains to show that the statement is true for a small system size, which is done by observing that for 
L = 1 both the left and the right hand side must be iQ □ 



4 Boundary passage probabilities 

When L is odd, all link patterns have an unpaired odd site. In the lattice this site belongs to an open 
path extending from — oo to oo. In this section we will calculate two probabilities: Pbithe probability 
that this infinite path passes through the first site at a given vertical position (Figure^); and Ph, the 
probability that this loop passes through the left boundary at a given vertical position (Figure ^p) . 

4.1 Definitions 

We first define to be the ground state eigenvector of the rotated system, given by 



a 



and related to l^*) by 

V5a(zi,...,Zi) = '(/'a(zL, . . . , Zl), 

where a and a are related by a rotation of tt. 



*Note that this proof is only vahd for odd L; for even L we must prove the statement for L = 2 as well. We omit this 
part of the proof as we are only interested in odd system sizes. 
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(4.1) 



Definition 4.1. The first site passage probability is given by 
where p acts between a link pattern and a rotated link pattern, 

giving 1 if the open path formed by these two link patterns goes through the first site, and if it does 
not. For example, the following configuration gives a weight of 1 (/> is depicted as two dots marking the 
first site): 

Definition 4.2. The boundary passage probability P^^^ is given by 




where p marks out the left boundary loop in the transfer matrix. It is depicted as 



(4.2) 























































































w 























Z; Z2 



ZL-7 ZZ. 



and like p, acts between and upward and a downward link pattern, multiplying each term in the transfer 
matrix by 1 if the infinite loop passes through the left boundary loop, and if it does not. For example, 
the following configuration is multiplied by 1: 





4 























As {P\ a) — 1, Vq;,/3, the denominator of P^^^'' and P^'^"* becomes 

(*| ^)=Y^ i^pi^a {P\ a) = ZLiZL, ...,Zi) Zl{zi, ...,Zl) = [Zl{zu- ■ ■ , ZL)f ■ 

a, 13 

Example: L = 3 

For L = 3 there are two link patterns, |() and ()|. Solving the qKZ equation for the components of the 

eigenvector gives 



tpQ\izi,Z2,Z3) 

and Z3 = V'K) + "^Ol is simply Xsizf, zl.zl). 



i>\{){zx,zi,z-i) = fc(z2,Zi), 

[9Z2/23] 1^ . 1 

k{l/Z2,Z3), 
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(3) 

The only combination of upward and downward link patterns which does not contribute to P^J is 
a — f3 ~ {)\. We thus find the first site passage probability by brute force as 

= ^2 {i^ioi^l, Z2, Z3)ijQi{z3, Z2, Zi) + tpQiizi, Z2, Z3)tpQi{z3, Z2, Zi) 
^3 

+tp\0{zi, 22, Z3)iIjiq{z3, Z2, Zl)) 

1 /'r I 4 I -i , r,/ 2 2 , 2 -2 I -2 2 I -2 -2 , 2 2 
= -^(5 + Zi+Zi +2(ZiZ2+ZiZ2 +Zi Z2+Z^ Z2 +Z1Z3 

I 2 -2 I -2 2 I -2 -2\ I 2 2 I 2 -2 , -2 2, -2 -2\ 
"T^l ^3 ~r ^1 Z3 -\- Z-^ Z3 j + Z2 ^3 + ^2 ^3 ' ^2 Z3 + Z2 Z^ j . 

4.2 Symmetries 

Proposition 4.1. P^^"^ is symmetric in Z2, ■ ■ ■ , zj^ and invariant under Zi — l/z^ for i > 2. 

Proof. This proof is similar to the proof that is symmetric, and uses the qKZ equation. Remem- 
bering that an i?-matrix acting between sites i ^ 1 and i + 1 commutes with p, we insert the identity 
Ri{zi+i/ Zi)Ri{zi/ Zi+i) into the definition for -PfJ^'': 



[L], ^ _ {'^L{zi,Zi+i)\Ri{zi+ilzi) p Ri{zi/ Zi+i)\'i L{zt, Zj+i)) 

, \. . . Zi, ZiJ^l . . . ) ~ 

ZL[Zi,Zi+iy 

(*L(^j+l,2:i)| Zi)) 



ZL{zi+i,ziy^ 

Pi "*(••• ^i+lj Zi . . .) . 



The invariance of pj^^"^ under Zi — > l/zi simply follows from the invariance of and under 

— >■ 1/zl, as well as the above symmetry. □ 

Proposition 4.2. -P(|^' is symmetric in all the Zi 's and invariant under Zi — )• 1/zi, Vi € {1, . . . , L}. 

Proof. The proof is similar to the previous proof, however it uses the fact that the interlacing condition 
(3.6 1 for the transfer matrix is also satisfied by p. □ 

4.3 Recursions 

Proposition 4.3. P^^'' satisfies the following recursions: 

Pl'^^lzlHiz.)^^- = Pl'''^\h,ZL), l<k<L. (4.3) 
Proof. The denominator of Pj^^'' has the recursion (3.12) 

L-2 

Zl{zI = q^zl_y^ = Y\ k{zL-i,Zif Zl-2(zi, . ■ ■,zl-2)'^, 

i=l 

and we will show that the numerator has the same recursion factor. 
From the recursion of the right eigenvector we have 

L-2 

P \zl = iqz^-ir ^ n (*i(4 = 1^zl_^)\p |*L-2(^1, • • • , ZL-2)) , 

i=l 

and since p commutes with (fL-i, we can consider {^l\(pl-i, which is the tt rotation of the vector 
(p\ \'^l{zl, ■ ■ ■ , zi)). Here, f\ is the bottom half of the TL operator ei, sending a link pattern of size L to 
one of size L ~ 2. We can thus obtain our desired result by calulating ei \^l{zl, ■ ■ ■ , zi)) and removing 
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the resulting link from site 1 to site 2. Here the gKZ equation comes in useful, as 



[zl-i/qzl] , [zl/qzl-i] 



[zl/zl-1 



zi=(92i.-i)^ 



[1/9] 



\^L{zL-l,qZL-l, ■ ■ -iZi)) 



L-2 



= |*L-2(^L-2, • • ■ , ^1)) X k{zL-l,Zi). 

Therefore, we have 

f\ \'^l{zl,---,Zi)) |4 = (q^^_,)2 = \^L-2{zL-2,---,Zi)) X fc(zi_i,Z,), 

which is the vr rotation of (^^-2! x Ili^i'^ ^C-zl-i, Zi), and thus, 

L-2 



L-2 



It follows that P, 



(i)l 



(L-2) 



The other relations in (4.3 1 follow from the invariance of Pi, under o Zj and under — > l/z^, for 



Proposition 4.4. -P^^'' satisfies the following recursions: 



^l4=(,..)±^ = Pb"^ ^'i^k,SL), l<k<L. 



_ p(i-2) 



□ 



(4.4) 



Proof. This proof is very similar to the previous one, but it also relies on the fact that p satisfies the 
same recursion as the transfer matrix (3.8), 



= qzi) o ip^ = (fi o p'^^ ^)(zi,Zj+i). 



□ 



4.4 Exact solution 

4.4.1 Inhomogeneous system 

Proposition 4.5. The explicit formula for the first site passage probability is 

p{L) _ XL-i{zl ■ ■ ■ , zl)xL+i{zl,zl z|, . . . , zl) 

' " xL{zi...,zir ■ ^'-'^ 

Proof. There are three steps to this proof. First, the degree of the proposed expression must be shown 
to agree with the definition. Second, enough recursions (or values of the polynomial at specified points) 
must be found to satisfy the degree. Thirdly, the proposed expression must be shown to be true for a 
small size example {L — 1), to initialise the recursion. 



The denominator of (4.1 ) can be easily shown to be Zf^, which has a degree of i — 1 in each variable 



zf , because L is always odd. The action of p can not raise the polynomial degree, so the numerator of 
pj^^^ must have at most the same degree as the denominator. The degree of the numerator of (4.5) is 
{L - 3)/2 + (L - l)/2 = L - 2 in each z^j*] 

*The definition of the transfer matrix implies that the components of the eigenvector, and thus the passage probabilities, 
are only functions of and do not depend on any odd powers of Zi. This fact is crucial to these proofs. 
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We view pj^ ^ as a polynomial in zf^ with coefficients in the other z^'s, and to satisfy the degree we 
need to know the value of the polynomial for at least 2L—1 values of zj^. The following recursions, which 



come from Proposition 4.3 will give the value of P^^"* at 4(i — 2) values of z£: 



Pt^U=(,z,)±- = ~'Hzfc,ZL), l<k<L. 



b lz£ = (<?Zfc) 

The proof that these recursions are satisfied by the proposed expression for -P(|^'' is straightforward, based 
on the known recursions for x (3.11 1. 

It thus remains to show that the expression is true for L — 1, which is easy to do as P^^"^ must trivially 
be 1, and = Xi = X2 = 1- □ 

Proposition 4.6. The explicit formula for the boundary passage probability is 

p{L) ^ -[g][w^] ^ x£+l(w^ z^ ■ . ■ , zl)xL+i{{q/wf,zl, ...,zl) 



Proof. This is a similar proof to the previous one, with the same three steps. 

The denominator of (4.2 1 is Z| multiplied by the denominator of p, which is the same as the denomi- 
nator of the transfer matrix; that is, Hi^i k{l/w, zi). Thus the degree of the denominator is L in each zf. 

Again because of the definition of p the numerator of P^'^'' will be at most the same as the denominator. 
The degree of the numerator of (4.6 1 is 2(L — l)/2 = L — 1 in each zf. 

We view P^^"* as a polynomial in z| with coefficients in the other Zi's and w, and to satisfy the degree 
we need to know the value of the polynomial for at least 2L + 1 values of z£ . As in the previous proof, 
we list here recursions which will give the value of P^^'^'' at 4(L — 1) values of z£. These recursions come 
from Proposition |4.4[ 

n''^^l4 = (g..)±^ = H''^'\h,ZL), l<k<L. 

The proof that these recursions are satisfied by the proposed expression for P^^^ is again based on the 
known recursions for x ( 3.11[ ), but also relies on properties of k(a, b) which imply that 

k{zk,w)k{zk,q/w) _ k{zk/q,w)k{zk/q,w/q) _ ^ 
k{l/w,Zk)k{l/w,qZk) k{l/w, Zk)k{l/w, Zk,q) 

It thus remains to show that the expression is true for L = 1, 

L -iq]i^^] ^ X2{w^,zl)X2{{q/w)^,zf) _ 



" k{l/w,z,) xMY k{l/w,z^y 

There are four configurations of p, three of which contribute to P^^^^ . The one which does not has a 
weight {[qzi / w][q/ ziw]){[qw / zi][qziw])~'^ , so 

p(l) [qzi/w] [q/ziw] ^ [qw/zi][qziw] - [qzi/w][q/ Ziw] ^ -[q] [w^] ^ 

[(7w/zi][gziw] k{l/w,zi) k{l/w,zi) 

4.4.2 Homogeneous limit 

In this section we will use the homogeneous limit zi = ... = z^ = 1, w'^ — —q, which causes the two 
orientations of the lattice faces to each have a probability of 1/2. 

Proposition 4.7. The expression for P^^^ in the homogeneous limit is 



piL) Av{L)Av{L + 2) 



where Av{L) and Ns(L) are the number of L x L vertically symmetric alternating sign matrices and 
cyclically symmetric self- complementary plane partitions of size L x L x L respectively, and have the 
explicit expressions 

'^^ (3» + 2)(6^ + 3)!(2» + l)! "f.^ (3» + 1)(6»)!(2»)! 

Ay{2m + 1) = [[ + 2)1(4^ + 3)! ' "^^^'"^^ = Ij (4z)!(4* + l)! ' 
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Proof. The result (4.7) is simply obtained from (4.5) and the homogeneous limit of the symplectic char- 
acters IT5]. 



X2™(l,...,l) = 3™('"-i)^y(2m + l), 
X2™-i(l,...,l) = 3('""')'^8(2m). 

Proposition 4.8. The expression for pj^^^ in the homogeneous limit is 



A{Lf 



4i7V8(L+l)My(i)2 ' 



□ 



(4.8) 



where A[L) is the number of L x L alternating sign matrices, A{L) — Y[i=o (l+iy. ■ 



Proof. For in (|4.6|, taking the homogeneous limit gives us 



pm ^ (-3) XL+i(-g,i,...,i)" 

(-4)i32(™-i)'Ar8(L + l)2' 
where we have set L = 2m — 1. To simplify the numerator, we use the relation |29j 

{4m — 2)/ —1 — 1 \ / \ / 

S\ , • ■ ■ ,M2m-l,'«2m-l) = X2m('7, Ul, • ■ ■ ,M2m-l) X X2m(-9,'"1, 



, U2m-l) , 



where s\ is the Schur function for the partition A. The homogeneous expression for s\ is 



4'^)(l,...,l) = 3^(^-i)/2^(L), 



and thus Xl+i{~<}, 1, ■ • • , 1) becomes 



XL+i{-q, 1, ■ • • , 1) 



3(2™-i)(™-i)^(i) 



Xl+i(9,1, ••■,!) 
We then use the recursion for XL+i{q, 1 • • ■ , 1) to get 

3(2m-l)(m-l)^(^) 



3(™-i)'A(L) 



fc(l,l)2'™-2;^2,n-2(l,...,l) Ay (L) 



from which the result (4.81 follows immediately. 



□ 



Proposition 4.9. In the limit L — ^ cx3, P^^^ and pj^^^ as given in (4.71 and (4.8) have the asymptotic 



4.4.3 Large-L limit 

Propositi 

behaviour 

where 



(4.9) 



9x 2-5/3r(i/3)r(5/6) 2i4/9^i/3G(l/3)4G(5/6)4 
r(l/6)r(2/3) ' G(l/2)4G(2/3)2 



Hence, the lattice probabilities Pt and Pf, follow the limiting behaviour (2.2 1 predicted by Schramm's 



formula at k = 6 with the correct exponent, but with different multiplicative constants. This is because 
Pf, and Pb correspond to fixing a position j and letting L tend to infinity, whereas the scaling limit would 
require a fixed ratio x = T^j / L. In the following section, we address the determination of P\cit for general 

3- 



5 Left-passage probabilities 

5.1 Definitions and properties 

Definition 5.1. For j g {1, . . . , i}, we denote by Xj{zi, . . . , zl) the probability that the path 7 passes 
through the j-th horizontal edge in a horizontal section of the type shown in Figure |2k. 
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Definition 5.2. For j G {0, . . . , L + 1}, we denote by Xj^w; zi, . . . , z^) the probability that the path 7 
passes through the j-ih horizontal edge in a horizontal section of the type shown in Figure [2]d. 

Definition 5.3. For j G {1, . . . , L + 1}, we denote by Yj{w; zi, . . . , zl) the probability that the path 7 
passes through the j-th vertical edge above a horizontal section of the type shown in Figure ^ 



Lemma 5.1. Using the convention that the vertices on a horizontal section like in Figure^Bp, are numbered 
{1/2, . . . , L + 1/2}, the probability that 7 passes to the left of {j + 1/2) is 



3 



1=1 

On a horizontal section like in Figure^), the probability that 7 passes to the left of {j + 1/2) is 

In particular, we have Xi — Pzji — Pb cind Xq = P1/2 = Yi ^ Ft- 

The power of (—1) in the above sums can be explained in the following way. The infinite loop, oriented 
as in Figure [2j can only pass in one direction through any given edge: For Xj, it passes upwards if j is odd 
and downwards if j is even; for Xj the opposite is true. When calculating the left-passage probabilities 
all configurations with the path passing downwards through the edge must be counted with a negative 
sign, thus we arrive at the above expressions. 

Lemma 5.2. The probabilities X , X and Y are related by the conservation property 

Vje{l,...,L}, Xj+X^=Y^+Y^+,, 
and the identities for special values of w 







J \W — Zj 








J \w—qzj 





Y, 



3 + i\w=qzj - ^3 



w=qZj 



Proof. The conservation property comes from considering the possible configurations on a face. Each X 
and Y has two terms corresponding to the two possible configurations, and for each X one of these terms 
matches one of those for one of the Ys. At the special values w — Zj and w — qzj, the face at position j 
is specialised to one of the configurations, and the identities follow immediately. □ 

Similarly to and Pb, the probabilities X, X and Y satisfy some symmetry and recursion relations. 

Proposition 5.3. The probabilities Xj{zi, . . . , zl) and Xj(w] zi, . . . , zl) are symmetric functions of 
{zi, . . . , Zj-i} and {zj+i, . . . , zl} separately. Xj{zi, . . . , zl) is invariant under zg — >■ l/zg for all I, and 
Xj{'w; 2i, . . . , 2i) is invariant under zg l/z£ for £ ^ j. The probability Yj{w; zi, . . . , z^) is a symmetric 
function of {zi, . . . , Zj_i} and {zj, . . . , z^} separately, and invariant under zg l/z£ for all i. 

Proof. The proof is similar to the ones for Proposition |4.1| and Proposition |4.2[ □ 

Proposition 5.4. The probabilities Xj{zi, . . . , z^) and Xj{w] zi, . . . , z^) satisfy the recursion relations 



X 



X 



(L) 
3 



X 



{L-2] 



X 
X 



(L) 
3 



{zi,Zi). 



X 



(L-2) 



{ze,ZL) ; 



forl<£< j, 
for j <i<L, 



whereas Yj(w; zi, . . . , z^) satisfies 



Y„ 



Y„ 



(L) 



rf-')(z,,zi). 



forl<£< j, 
for j < i < L. 



'We could also define Yj in a similar fashion, but this is simply related to Yj by the transformation l/^fc, Vfc. 
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Proof. The proof is similar to the ones for Proposition |4.3| and Proposition |4.4[ 



□ 



The above relations satisfied by Xj, Xj and Yj are, in principle, sufRcent to determine completely 
these quantities. However, they turn out to be particularly difficult to solve in practice, and we resort to 
different tools to describe them. 



5.2 Numerical study 

For a given choice of the Zj 's, one can find numerically the components ipa by the power method, and 
use these to evaluate Xj, Xj and Yj. The left-passage probabilities are then obtained from Lemma 



5.1 



For small enough system sizes, we find the left-passage probabilities Pj^i/2 as rational fractions (see 
Table [l]). 



L 


Zl 


Pl'^ X Zl 


3 


2 


0, 3, 1, 4 


5 


11 


0, 78, 22, 99, 43, 121 


7 


170 


0, 16796, 4484, 21093, 7807, 24416, 12104, 28900 


9 


7429 


0, 29641710, 7721790, 37074705, 12859293 . . . 


11 


920460 


0, 426943865250, 109785565350, 532943651700, 178807268772, 605036201854 . . . 



Table 1: Left-passage probabilities Pj^^ for strips of width L < 11. Data obtained by numerical 
diagonalisation of the transfer matrix. 



We now turn to the convergence of Pj+1/2 to Schramm's formula (2.1 ). From Lemma 5.1 we see that 
Pj+1/2 is the sum of an alternating sequence, and has oscillations of wavelength Sj~l. This phenomenon 
appears clearly in Figure [4j For this reason, we define the smooth and oscillatory parts as 



P. 



Pj I {Pj-i/2 - Pj+1/2) = li-iy-'x, 



In Figure [4] we see that Pj is very close to Schramm's formula for L — 21. In Figure [5] we compare the 




Figure 4: The left-passage probability Pj+1/2 and its smooth part Pj for L = 21, compared to 
Schramm's formula (2.1|. 



data for Pj at various system sizes, and observe very good convergence to Schramm's formula. Finite-size 
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Figure 5: Smooth part of the left-passage probability, compared to Schramm's formula (2.1) 



effects are more important near the boundaries, but as we already noted in Section |4j the scahng of P1/2 
with L is the one predicted by Schramm's formula. 

Finally, in Figure [6j we plot the oscillatory part Pj. This quantity is a lattice effect, and is not 
predicted directly by Schramm's formula. 



5.3 Left-passage probability in the FK model 
5.3.1 Mapping to the six-vertex model 

Throughout this section, we remove the restriction on g, and we use the algebraic Bethe ansatz notations 
Moreover, we introduce for convenience 77' = + i7r, so that the loop weight reads 



-2 cosh r; = 2 cosh 77' . 



Following 1^, we distribute the loop weight locally by orienting the loops: to each 7r/2 left (resp. right) 
turn of a loop, we associate a phase factor (resp. in the Boltzmann weights. Forgetting 

about the loop connectivities results in a six- vertex (6V) model (see Figure [t]). Using the Boltzmann 
weights defined by the i?-matrix in Definition 3.1 (with z — I wlog for the moment), the 6V weights can 
be rescaled to 

sinh(77 -|- u) 
sinh u 

e+^'/s+^sinhry ^^'^^ 
Moreover, on the boundary, the loops undergo a half-turn, and hence the boundary weights must be 



W3 



UJ2 

We 



-7772 



Hence the 6V model resulting from this mapping is described by the matrices 



(5.2) 



R{u) 



/ sinh(77 




u) 



V 





sinhu 
6^+" sinh r] 






sinh 77 

sinhw 









sinh(77 - 



u) J 



K± 
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and the corresponding transfer matrix is denoted by tgy. These matrices are related to the standard 6V 
ones by the "gauge transformation" : 



where K+{u) = 2e«+if (71 + ?7, K^{u) = 2e«- ^„), and 



(5.3) 



R{u) 



/ sinh(?7 + u) 

sinhit sinh?7 

sinh77 sinhw 

\ sinh(r; + u) / 



sinh(^ + u) 

sinh(^ — u) 



with the values of the boundary parameters: ^± = =Foo. It is customary to shift the spectral parameters 
to define the monodromy matrices 1251 : 



:= A - 77/2 , V, := ~ 77/2 , 
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and we thus write: 



r(A) 

f(A) 
i6v(A) 



= i?oL(A-eL)---^oi(A-ei) 
= Tro \K+{X)T{X)K_{X)f{\) . 



(5.4) 



One can easily show that the transfer matrices before and after the gauge change (5.3 1 are simply related 
by a similarity transformation 



^6V — ^ ^ ^6V G , 



where 



3 = 1 



(5.5) 



If we specialise to a homogeneous system where all the ^j's are set to 77/2, the very anisotropic limit 
A 77/2 yields the open XXZ Hamiltonian 



dlogteviX) 



dX 



L-1 



E [^j + ^J'^J+i + coshry c^|cr|+i] - sinhT] (df - fjf,) 



(5.6) 



Note that in the critical regime (ry S iM), the boundary terms are imaginary. In the remainder of this 
section, we shall restrict ourselves to the homogeneous system described by 'Hxxz, but our results can 
be readily generalised to an arbitrary choice of the fj 's. 



5.3.2 The left-passage probability as an XXZ correlation function 
Proposition 5.5. In the critical regime 77 G iM, the following identity holds: 



P. 



j+1/2 



5: Re 



i=i 



(*o 1*0) 



(5.7) 



where and o.tr the left and right eigenvectors ofHxxz (5.6| associated to the lowest energy. 

Proof. We have used above the mapping of TL loop configurations to 6V arrow configurations, through 
the orientation of loops. Consider the intermediate model, i.e., the oriented TL (OTL) loops. In this 
model, a configuration C on the whole lattice has a Boltzmann weight (at the isotropic point A — — 77r/2) 

j^j'^j ^ ^^#anti— clockwise loops(C) ^ ^ ^clockwise loops(C) 

and the partition function is equal to the one of the original TL model 

oriented config. C unoricntcd config. C 

We now consider a horizontal section of the strip, of the type shown in Figure [2^. For a given oriented 
loop configuration C, we denote <7|(C) £ {1, —1} the orientation of the arrow across the j-th horizontal 
edge. There are three possibilities for the loop passing through this edge (see Figure [s]) : 

• The loop is the open path 7. Then o-|(C) ~ (—1)-'^^. 

• The loop encloses the marked point on the left of j. Then crj{C) = +1 iff the loop is oriented 
anti-clockwise. 

• The loop encloses the marked point on the right of j. Then <Jj{C) — +1 iff the loop is oriented 
clockwise. 

We denote by Aj,A'j,Aj the corresponding subsets of oriented loop configurations. The expectation 
value of crj in the OTL model then reads 



1 

Z 



(-1)'"' E E - E 



ceAj 



^ceA', ceA'fj 



W\C] 



i~iy-^ Xj + tanhr/ (P[C G MA - P[C G M']) 



(5.8) 
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Figure 8: From left to right: loop configurations belonging to the subsets Aj, A'j, and Aj, for 



2 and L 



Finally, from the mapping described above, we have 

(*o|*o) 



□ 



6 Perspectives 

In the percolation model, we have obtained an exact expression for the boundary passage probabilities 
Pf, and Pb, which are lattice analogs of a boundary observable in SLEg. In the FK cluster model with 
generic Q, we have related the left-passage probability to the magnetisation in a solvable open XXZ spin 
chain. 

Our results have many possible developments. First, within the qKZ approach, we hope to exploit 
the symmetry and recursion relations for the probabilities Xj, Xj and Yj to find their explicit expres- 
sion. This certainly involves a deeper understanding of the properties of symplectic characters and Schur 
functions [S] . Second, with the algebraic Bethe ansatz [T71 [TS] , it seems possible to calculate the magneti- 
sation ((t|) in the open XXZ chain in a closed form, at least for j close enough to one of the boundaries. 
The advantages of this method is that it is valid for any value of the deformation parameter q, and, in 
cases where a closed form cannot be achieved, it still produces determinant forms which can be evaluated 
numerically for very large system sizes (L ^ 1000 sites). Finally, we note that the probabilities Xj, Xj 
and Yj are very similar to the discretely holomorphic parafermions found for the TL loop model [161 122) , 
which are the starting point in the proof of conformal invariance for the Ising model [HI [IZ] ■ Thus the 
study of these objects on a general domain may bring progress in extending this proof to the FK model 
with generic Q. 
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